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Abstract 

Complete eigenfunctions of linearized integrable equations expanded around an arbitrary solution 
are obtained for the Ablowitz-Kaup-Newell-Segur (AKNS) hierarchy and the Korteweg-de Vries 
(KdV) hierarchy. It is shown that the linearization operators and the recursion operator which 
generates the hierarchy are commutable. Consequently, eigenfunctions of the linearization oper- 
ators are precisely squared eigenfunctions of the associated eigenvalue problem. Similar results 
are obtained for the adjoint linearization operators as well. These results make a simple connec- 
tion between the direct soliton/multi-soliton perturbation theory and the inverse-scattering based 
perturbation theory for these hierarchy equations. 



1 Introduction 



Integrable equations are nonlinear evolution equations which can be solved exactly by the in- 
verse scattering method. Over the past few decades, it has been discovered that many physi- 
cally important equations such as the Korteweg-de Vries (KdV), nonlinear Schrodinger (NLS) and 
sine-Gordon equations are integrable (see [1] and the references therein). Linearization of an inte- 
grable equation around its solution arises in many important applications, most notably in a direct 
soliton/multi-soliton perturbation theory. In such situations, eigenfunctions of linearization opera- 
tors and their completeness are the key questions. For linearization around single-soliton solutions, 
these complete eigenfunctions have been obtained for a large class of integrable equations such as the 
KdV hierarchy, NLS hierarchy, modified-KdV hierarchy, sine-Gordon, and Benjamin-Ono equations 
[2, 3, 4, 5, 6, 7, 8, 9]. It has been found that these eigenfunctions are related to squared eigenfunc- 
tions of the associated eigenvalue problem (except for the Benjamin-Ono equation). However, for 
linearization around a general solution such as a multi-soliton solution, complete eigenfunctions are 
known for much less integrable equations [2, 8]. Some general ideas have been proposed to deter- 
mine these eigenfunctions though. One idea by Keener and McLaughlin [2] is that eigenfunctions 
of a linearization operator expanded around an arbitrary solution are the variations of the solution 
with respect to each parameter in the scattering data. Another idea by Herman [6] is to utilize 
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the Lax pair of the integrable equation and find special combinations of squared eigenfunctions 
of the associated eigenvalue problem, so that these combinations satisfy the linearized equation of 
the evolution equation. But in both approaches, each equation has to be treated separately. In 
addition, for each equation, much work is needed to find eigenfunctions of the linearization oper- 
ator, or relate them to squared eigenfunctions of the associated eigenvalue problem. The idea by 
Yang [9], however, is free of these problems. This idea is to show that linearization operators of 
a hierarchy and the recursion operator which generates the hierarchy are commutable, thus they 
share the same set of eigenfunctions. Furthermore, these eigenfunctions are simply squared eigen- 
functions of the associate eigenvalue problem. Compared to the other two approaches, this method 
explicitly gives the eigenfunctions of linearization operators in the simplest way. In addition, it 
treats an entire hierarchy all at once. In [9], this idea was applied only to linearizations of the KdV, 
NLS and modified-KdV hierarchies around single-soliton solutions. In that special case, our results 
went beyond commutability of the operators. We also showed that linearization operators of the 
hierarchy equations could be factored into the recursion operator of the hierarchy and the lineariza- 
tion operator of the lowest-order equation in the hierarchy. Commutability between linearization 
operators and the recursion operator is a simple consequence of this factorization representation 
for the linearization operators. 

In this article, we extend the results of Yang [9] to the Ablowitz-Kaup-Newell-Segur (AKNS) hi- 
erarchy and KdV hierarchy linearized around an arbitrary solution. The AKNS hierarchy is the 
family of integrable equations associated with the Zakharov-Shabat eigenvalue problem [10, 11], 
and the KdV hierarchy is associated with the Schrodinger eigenvalue problem [11, 12]. In this 
general case, we can still show that linearization operators are commutable with the recursion 
operator of the hierarchy (the factorization result of linearization operators for single-soliton so- 
lutions does not extend to this general case though). This commutability allows us to establish 
that complete eigenfunctions of linearization operators in the AKNS or KdV hierarchy are simply 
squared eigenfunctions of the Zakharov-Shabat or Schrodinger operator. Similar results can be 
obtained for adjoint linearization operators as well. In a direct soliton/multi-soliton perturbation 
theory, these squared eigenfunctions will then serve as the expansion basis for perturbation so- 
lutions. Interestingly, these same squared eigenfunctions were also used to expand perturbation 
solutions in the inverse-scattering based perturbation theory [13, 14, 15]. Thus, our results in this 
article indicate that, at a deeper level, the direct soliton/multi-soliton perturbation theory and the 
inverse-scattering based perturbation theory are really equivalent. 



2 Eigenfunctions of linearization operators for the AKNS hierar- 
chy 

The AKNS hierarchy associated with the Zakharov-Shabat eigenvalue problem is [11]: 



n 



+ iu(2Lt) 



= 0, 



(2.1) 



where the recursion operator L+ is 



2i 
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l + 2 q f* oc dyr 



(2.2) 
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and co(k) is the dispersion relation of the linear equation in the r-component. In this section, we 
require that co(k) is an entire function. When uj(k) = k 2 and q = —r*, Eq. (2.1) reduces to the NLS 
equation; when uj(k) = k 3 and q = — r, the modified-KdV equation results. The adjoint operator 



of L+ is 



L z = 



2i 



- m ~ 2 i iT d v r - 2 <iir d y<i 

2r dyr m + 2r XT d V1 



(2.3) 



Now suppose [ro(x, t), qo(x, t)] T is an arbitrary solution of the evolution equation (2.1). To avoid 
dealing with divergent integrals in the following analysis, we require that this solution vanish as 
\x\ goes to infinity. But generalization to non- vanishing solutions is also possible by appropriately 
defining divergent integrals, as we did in [9]. Next we linearize Eq. (2.1) around this arbitrary 
solution. For this purpose, we write 



r (x, t) + f(x, t), q = q (x, t) - q(x, t), 



(2.4) 



where f,q <S 1. Note that we deliberately introduced opposite signs in front of r and g's pertur- 
bations. This is important for obtaining the commutability relations which we will present later in 
this section. When Eq. (2.4) is substituted into Eq. (2.1), linearization of Eq. (2.1) is: 



L 



0, 



(2.5) 



where L is the linearization operator. We denote the adjoint operator of L as L A . We also denote 
Lq and Lq as the recursion operators Lf and L z with r and q replaced by the solutions ro(x,t) 
and qo(x,t). The primary objective of this section is to show that operators L and Lq commute, 
and L and Lq commute, i.e., 

LL+ = Lq L, (2.6) 



and 



L A Lq — LqL A . 



(2.7) 



We will prove relation (2.6) first. Relation (2.7) will then follow naturally. 



Without loss of generality, we assume that the dispersion relation u>(k) is a power function, uj(k) = 
k n , where n is a non-negative integer. The reason is that any entire function of uj(k) can be 
expanded into a power series. Linearization of operator Lf around the solution (ro,<7o) is 



+ 0(f\?q,q 2 ), 



(2.8) 



where the operator T is defined as 



r 
Q 



-2f JZoo dyqo + 2r /f ^ dyq 2? /f ^ dyr + 2r /f ^ dyr 
2q /f oo dyqo + 2q /f dyq -2q /f dyr + 2q /f dyr 



(2.£ 



Then, for power functions of co(k), the linearization operator L is simply: 



L 



+<(2L+r 



+ ^(2L+)^ 



fc=i 



(2L + ) 




(2.10) 
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Denoting 



Pa 
Qn 



ro 



then the evolution of (ro,qo) becomes 



-qo 



Pn 
Qn 



(2.11) 



(2.12) 



and the functions (P n ,Qn) satisfy the recursion relation 



Pn+l 
Qn+l 



2L r 



P 
1 n 

Qn 



(2.13) 



To establish the commutability of operators L and L^, we examine the function 

H n ^2z(L+L-LL+)(^ 
Simple calculations show that H n has the following expression: 
7~L n = 



(2.14) 



2-Pn f-ooiqor - r q)dy - 2r ^{Qn? + P n q)dy 
-2Qn I-ooiqor - r q)dy - 2q J^iQnf + P n q)dy 



-(2L+ 



\n+l 



+ (2^ + 



-oo \ 
r 

-oo v 

1 

-1 



r 

q 



+ ELi(2^ + ) fc - 1 teat? 



T 



2iU 



r 

q 



(2L+ 

(2L+r~ k 



qo 



(2.15) 



Here we have replaced the time derivatives rot and —qot in Ti n by P n and Q n in view of Eq. 
(2.12). It is important to realize that the above 7i n expression (2.15) is now purely algebraic, and 
is independent of the evolution equation (2.12). Below we will use algebraic manipulations and the 
induction method to prove that TL n is zero for all n > 0. 

When n = or 1, one can verify directly that 7i n is indeed zero. Now we assume that 7i n = 
for some n > 0. Then we try to show that H n +i = 0. For this purpose, we calculate the quantity 
Hn+i — 2LqHu. It turns out that most of the summation terms in 7Y n +i and 2Lq"H„ cancel each 
other out. The terms left over are 



W n +i — 2L§H n 



2P n+ i Xf ^ (q r - r q)dy - 2r f*^ {Q n+1 f + P n+1 q)dy \ 
-2Q n+1 jl^qof - r q)dy - 2q JZ^Qn+if + P n +iq)dy J 

J-ooiQnf + P n q)dy \ 
2 qo J-ooiQnf + P n q)dy J 



-9T+ ( 2P n J-oo(qor - r q)dy - 2r , 
^° { -2Q n f* OD (q f-r q)dy-2q( 



2L+F 



2iL+ 



p 

1 n 
Qn 



(2.16) 



When the recursion relation (2.13) for P n +i and Q n +i is substituted into the above expression, 
algebraic simplifications immediately reveal that 



'Hn+i — ZLqHu — 0. 



(2.17) 
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Since H n is zero by assumption, it then follows that H n+ i = 0. Thus H n = for all n > 0, which 
means that L and Lq are commutable. For a general entire function of the dispersion relation u(k), 
this result still holds, as an entire function can be expanded into a power series. 

The proof for the commutability of L A and Lq is trivial once the commutability of L and Lq has 
been established. The adjoint operator of LLq is LqL a , and the adjoint of LqL is L a Lq. Since 
LLq = LqL, their adjoints are certainly the same, i.e., L L A = L A L . Thus L A and Lq are also 
commutable. 

An important consequence of the commutability relations (2.6) and (2.7) is that L (L A ) and Lq 
{Lq) share the same set of eigenfunctions. To see how this comes about, let us assume that *ff(x, t, Q 
is a continuous eigenfunction of Lq with real eigenvalue £, he., 

L+* = C*. (2.18) 

Under the condition that [ro(x, t),qo(x, t)] vanishes as \x\ goes to infinity, we can impose the bound- 
ary condition for ^ as 



-e~ 



C) — ( „_2ifx > x^-oc. (2.19) 



Since L and Lq are commutable, we have 

LqL^ = (Lff. (2.20) 

Thus L^ff is also an eigenfunction of Lq with eigenvalue Q. As x goes to infinity, the linearization 
operator L becomes 

/ d t + iuj(-id x ) \ || , . 

^ — M n « ■ / ■ q \ I 5 x -> oo. (2.21) 

^ d t -iuj(id x ) J 

Consequently, the boundary condition for L^ff can be obtained from Eqs. (2.19) and (2.21) as 

LV(x, t, C) — -iw(2C) ( _ e ° 2iCx ) , x ^ -oo, (2.22) 

which is proportional to the boundary condition (2.19) of eigenfunction iff. Then it becomes clear 
that L^f and * are the same eigenfunction of operator Lq with eigenvalue £ (i.e., they are linearly 
dependent). In view of their boundary conditions, we see that 

L^ff = -iw(2C)tf, (2.23) 

i.e., %f(x,t,Q is also a continuous eigenfunction of operator L with eigenvalue —iu(2Q. 

For the same real eigenvalue £, Lq has another linearly independent eigenfunction iff with boundary 
condition 

/ 2i(x \ 

*(M,C) — ( J, x^-oo. (2.24) 

Similar analysis shows that iff is also a continuous eigenfunction of L, but with eigenvalue iu(2Q, 
i.e., 

L^f = iw(2C)^. (2.25) 
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For the discrete eigenfunctions and generalized eigenfunctions of Lq~, same analysis indicates that 
they are also discrete eigenfunctions and generalized eigenfunctions of L. Thus Lq and L indeed 
share the same set of eigenfunctions. Naturally, the same statement applies to Lq and La as well. 

What exactly are the sets of eigenfunctions for L and L A 1 Are these sets complete? In view of 
our results above, we only need to find the answers for operators Lq and Lq. The eigenfunctions 
for Lq" and Lq and their closure have been known for over twenty years from the celebrated work 
by Ablowitz, Kaup, Newell and Segur [11] and by Kaup [16]. The results can be summarized as 
follows. 



Consider the Zakharov-Shabat eigenvalue problem with potential [qo(x,t),ro(x,t)]: 

vix + i£vi = q (x,t)v 2 , 
V2x ~ Kv2 = r (x,t)vi, 

and define Jost functions for real ( as 



ip(x,t,Q -- 
tp(x,t,() = 

4>(x,t,C) = 



■01 







1p2 




1 






' 1 


1p2 







01 " 




" 1 


02 







01 " 







02 




-1 



x — > oo, 



X 



oo, 



x — > — oo, 



x — > — oo. 



(f){x,t,() = 
The right and left solutions are related by 

<t>(x, t, o = a (t, cmx, t, o + b(t, ov(x, t, o, 

0(x, t, () = -a(t, CMx, t, C) + b(t, CWx, t, C), 

where 

a(t,Qa(t,Q+Kt,OKt,Q = l 
from Wronskian relations. With Eq. (2.34), the inverse of Eqs. (2.32) and (2.33) is 

■0 = — a0 + 60, 

-0 = a0 + 60, 



2.26) 
2.27) 

2.28) 

2.29) 
2.30) 
2.31) 

2.32) 
2.33) 

2.34) 

2.35) 
2.36) 



where we have suppressed the dependent variables x, t and £. In addition to the continuous spectrum 
(C real), Eqs. (2.26) and (2.27) may also possess discrete eigenvalues (bound states) in the upper 
and the lower half (-plane. In the upper half plane, these occur whenever a(t, Q = 0, and we 
designate them by £fc, k = 1,2,...,N, where N is the total number of bound states in the upper 
half ("-plane. At ( = (fcj and ip become linearly dependent and 



(x, t, Cfc) = b(t, (Mx, t, C fc ), k = l,2,...,N. 



(2.37) 
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In the lower half (-plane, bound states correspond to zeros of a(t, Q which we designate by £ k , k 
1,2,..., TV. At ( = Ck, 



(x, t, C fc ) = b(t, (kW(x, t, Cfe), k = l,2,...,N. 



(2.38) 



It is important to note that when [qo(x,t),ro(x,t)] is a solution of the AKNS hierarchy (2.1), the 
discrete eigenvalues ( k and Ck are independent of time t. 

With the above notations, the eigenfunctions and generalized eigenfunctions of operators Lq and Lq 
are simply squared eigenstates of the Zakharov-Shabat system (2.26) and (2.27) [11]. Specifically, 
the set of eigenfunctions and generalized eigenfunctions for Lq is 





' ft 










d 




! 






[ ft ' 


, ( real; 


[ ft ' 


\ ft 1 




. ~ft . 


C 


. -ft . 


C 


. -ft . 


. -ft . 



,1 < k < N; 



ft ' 


d 


ft 


. -ft . 




. -ft . 



and the set of such eigenfunctions for Lq is 



1 < k < N 
(2.39) 



ft 



ft 
ft 



, £ real; 



ft 
ft 



d_ 

Cft d ( 



ft 
ft 



,1 < k < N: 



\ ft ' 




d 


\ 1 


. ft . 






. ft . 



,1< k < N 



(2.40) 

It has been shown by Kaup [16] that each of these two sets is complete. The orthogonality and inner 
products of functions in these sets have also been obtained there. In view of these facts, we then 
conclude that the sets (2.39) and (2.40) are also the complete sets of eigenfunctions and generalized 
eigenfunctions for linearization operators L and L A respectively. What about the corresponding 
eigenvalues? The eigenvalues are actually quite easy to obtain from the asymptotic behaviors of 
these eigenfunctions. For operator L, the results are: 



and 



L^ 



L 



L 



L 



L 



-ft 



Ck 




= -iw(2C) 

= iu(2() 
-iu(2( k 



Ck 



= -iuj(2( k ) — 



-ft 



ft 
-ft 



ft 



Ck 



Ck 

2ic/(2Cfc) 



Ck 



C real; 



C real; 



1< k < N: 



1 < k < N; 



Ck 



1< k < N: 



(2.41) 

(2.42) 
(2.43) 
(2.44) 
(2.45) 



Ck 



- d 
iu{2( k ) — 



+ 2iJ{2Q k ) 



Ck 



Ck 



1 < k < N. 



(2.46) 



The results for L A are: 



L J 



-ilo(2Q 



ft 
ft 



C, real; 



(2.47) 
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and 



L dC 



L J 



d_ 



M2() 
-iu(2( k ] 
iu(2( k ) 



ft 



C real; 
1 < k < N; 
1 < k < N; 



ft 
ft 



ft 
ft 



d_ 



ft 
ft 



■ 2iw'(2C fc ) 
+ 2ia/(2Cfc) 



ft 
ft 



Cfc 



ft 
¥2 



1 < k < N: 



1< k < N. 



(2.48) 
(2.49) 
(2.50) 
(2.51) 

(2.52) 



Cfc 



Lastly, we note that in the development of a direct soliton/multi-soliton perturbation theory, it 
is often convenient to use the derivatives of soliton/multi-soliton solutions with respect to soliton 
parameters as discrete eigenfunctions and generalized eigenfunctions of the linearization operator 
[2, 5, 9]. These derivative states span the same linear space as the discrete eigenfunctions in the 
set (2.39) do. Thus use of either discrete set is sufficient. 



3 Eigenfunctions of linearization operators for the KdV hierarchy 



For the KdV hierarchy, similar results hold. The analysis is simpler though as we only have a scaler 
equation to consider. This hierarchy can be written as [11]: 

q t + C(ALf)q x = 0, (3.1) 

where q(x, t) is a real function, C(k 2 ) is the phase velocity of the linear equation, and the recursion 
operator Lf is: 

Id 2 1 f°° 

i ' + = -4a?-" + 2*i *■ (3 ' 2) 

Here the subscript "s" in Lf refers to "Schrodinger" , as the associated eigenvalue problem for the 
KdV hierarchy (3.1) is the Schrodinger equation [11, 12]. In this section, we require the phase 
velocity function C(z) to be entire. The adjoint operator of Lf is: 

Id 2 1 f x 

Ls = -4d^- q + 2Lj yq y (3 - 3) 



Suppose qo(x,t) is an arbitrary solution of the evolution equation (3.1). To linearize Eq. (3.1) 
around this solution, we write 

q = q (x,t) + q(x,t), (3.4) 

where g « 1 When Eq. (3.4) is substituted into the evolution equation (3.1) and higher order 
terms in q neglected, the linearized equation is 

L kh q = 0, (3.5) 
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where L k h is the linearization operator. Here the subscript "kh" is the abbreviation of the KdV 
hierarchy. The adjoint operator of L k h will be denoted as L kh . We will also denote Lf and L s q as 
the operators Lf and L s with q(x, t) replaced by the solution qo(x, t). The objective of this section 
is to show that Lkh and Lf are commutable, and L kh and L s q are commutable, i.e., 

L kh L+ = L+ L kh , (3.6) 

and 

(3.7) 

These results are analogous to Eqs. (2.6) and (2.7) for the AKNS hierarchy. 

Without loss of generality, we will just prove relations (3.6) and (3.7) for power functions of the 
phase velocity function, C(z) = z n , where n is a non-negative integer. For this power function, it 
is easy to check that the linearization operator L k h is: 

L kh q = qt + (4L+) n q x + f^^)*" 1 ^^^ (3.8) 

k=i 

where 

F kh [q] = -4q + 2q x [°° dy, (3.9) 

Jx 

and qo(x,t) is a solution of the evolution equation (3.1). Denoting 

W n = -(4L+) n q 0x , (3.10) 

then got is simply 

got = W n , (3.11) 
where functions W n satisfy the recursion relation 

W n+1 = 4L+W n . (3.12) 

To show that L k h and L^}~ are commutable, we calculate the quantity 

J n = -4(L+L kh - L kh Lf )q, (3.13) 

which has the expression 

Jn = 2W nx f™ qdy - 4W n q - (4L+) n [4L+q x - 4(L+q) x 

~ ELl(4L s + ) fe_1 {^f T kh [q} - J^[AL+g\} (4Lf r- k q 0x . (3.14) 
Now we show that J n = for all n > by the induction method. 

When n = or 1, trivial calculations show that J n is indeed zero. Now assume that J n = for 
some n > 0. Notice that 

J n+1 - 4L+J n = 2W n+1:X qdy - 4W n+1 q - 4L+ [2W nx qdy - 4W n q] 

+ {4L+F kh [q} - F kh [4L+q]}w n . (3.15) 
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Substituting the recursion relation (3.12) for W n+ \ into the above equation (3.15) and carrying out 
some algebraic simplifications including integration by parts, we find that 



Jn+l — ^L^Jn 



0. 



(3.16) 



Since J n = by assumption, we see that J n +\ = 0. This induction procedure proves that J n = 
for all n > 0. Thus, L^h and Lf are commutable for any power function of the phase velocity 
C{z). The commutability for general entire functions of C(z) follows from the fact that an entire 
function can be expanded into a power series. 

Now that and L^ are commutable. Taking the adjoint of the commutability relation (3.6), we 
find that and L s q are also commutable. 

Commutability of L^h (L^ h ) and Lf (L s q) implies that these operators share the same set of 
eigenfunctions and generalized eigenfunctions. The eigenfunctions of L^ and L s q are well known 
[11, 15]. They are simply squared eigenfunctions of the Schrodinger operator with potential qo(x, t). 
Specifically, consider the Schrodinger equation 



+ 



( 2 +qo(x,t) 



0. 



(3.17) 



Using conventional notation, we define the eigenstates ip(x,t,() an d 4>(x,t,() of Eq. (3.17) as 

-oo; 



e s , x — > oo; 

a(t,()e^ x -b{t,-Qe~< x , x 



(3.18) 



and 



<f>(x,t,() 



e %c - x , x — > — oo; 

a(t,C)e-^ x + b(t,C)e^ x , x -> oo. 



(3.19) 



In addition to the above continuous spectrum (real £), Eq. (3.17) may also possess discrete eigen- 
values in the upper half ("-plane (on the imaginary axis for real potential qo) where a(t, Ck) = 0,k = 
1,2, . . . , N . Note that if qo(x, t) is a solution of the KdV hierarchy (3.1), then these discrete eigen- 
values Cfc are independent of time t [11, 12]. With the above notations, the set of eigenfunctions 
and generalized eigenfunctions for the operator Lf is 



d^ 1 



dx 



, £ real; 



dx 



dV 2 



dxd( 



,1 < k < N 



(3.20) 



and the set of such eigenfunctions for L s q is 

, C real; <p 2 



,1 < k < N 



(3.21) 



Commutability of Lf Q (L s q) and Lkh (-^kh) shows that the sets (3.20) and (3.21) are also eigenfunc- 
tions and generalized eigenfunctions for the linearization operators L]<h an d respectively. In 
addition, we can readily show that the eigenvalue relations are 



dj> 2 
dx 



= 2iCC(4C 2 



dj) 2 
dx 



( real; 



(3.22) 
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L 



kh 



d^ 1 



dx 



di> 2 



dx 



d 2 ^ 2 



dxd( 



2 , av 2 



dxd( 



+ 



2iC(4( 2 ) + m( 2 k C'(4( 2 k 



1 < k < N; 
d^ 2 



and 



L 



Cfc 

A A? 
kh 9 



dx 



Cfc 



£kh0 2 



Cfc 



. = 2<C(4C 2 
2< fe C(4C|) , 



C' 



C real; 
1 < k < N; 



^kh 



= 2i( k C(4C 2 k ) 



[2iC(4$) + lGi$C'(4&j 



Cfc 



(3.23) 

1 < k < N; (3.24) 

(3.25) 
(3.26) 



1< k < N. 



(3.27) 



The completeness of the two sets (3.20) and (3.21) and their inner products have been derived in 
[15, 17]. Thus these sets can be used to expand the perturbation solutions in a direct soliton/multi- 
soliton perturbation theory [6, 7, 9]. 



4 Concluding remarks 



In this article, we have studied the linearization operators of the AKNS hierarchy and KdV hierarchy 
equations expanded around an arbitrary solution. We have found that these linearization operators 
and the recursion operator which generates the hierarchy are commutable. This commutability 
relation immediately reveals that linearization operators and the recursion operator share the same 
set of eigenfunctions, and these eigenfunctions are simply squared eigenfunctions of the Zakharov- 
Shabat or Schrodinger equations. Compared to the other methods for determining eigenfunctions of 
the linearization operators [2, 6], our method is simple, and it gives the eigenfunctions for the entire 
AKNS and KdV hierarchies all at once. In addition, our result makes a clear connection between 
the direct soliton/multi-soliton perturbation theory and the inverse-scattering based perturbation 
theory, as perturbation solutions in both theories are expanded onto the same complete set of 
eigenfunctions. With the eigenfunctions of linearization operators now available, one can proceed 
to develop a direct soliton/multi-soliton perturbation theory for the AKNS and KdV hierarchies, 
which should reproduce the results of [14, 15] obtained by the inverse-scattering based perturbation 
method. This problem falls outside the scope of the present article. Another interesting question is 
whether the idea of this paper can be extended to derive eigenfunctions of linearization operators 
for other integrable equations. We believe the answer is positive. For instance, our results should 
be extendible to the integrable N-w&ve equations whose recursion operator and its eigenfunctions 
have been available [18, 19]. This question will be left for future studies. 
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